On 
O 
00 



o 
o 



>< 



cd 



PACS numbers: 98.65,98.80 
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The halo-galaxy lensing correlation function or the average tangential shear profile over sampled 
£SJ \ halos is a very powerful means of measuring the halo masses, the mass profile, and the halo-mass 

correlation function of very large separations in the linear regime. We reformulate the halo- galaxy 
lensing correlation in harmonic space. We find that, counter-intuitively, errors in the conventionally 
used flat-sky approximation remain at a % level even at very small angles. The errors increase at 
r- — , larger angles and for lensing halos at lower redshifts: the effect is at a few % level at the baryonic 

r\l ■ acoustic oscillation scales for lensing halos of z ~ 0.2, and comparable with the effect of primordial 

non-Gaussianity with /nl ~ 10 at large separations. Our results allow to readily estimate/correct 
for the full-sky effect on a high-precision measurement of the average shear profile available from 
\^ • upcoming wide-area lensing surveys. 
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Dark matter halos hosting galaxies and galaxy clusters emerge as a result of gravitational amplification of tiny 
primordial fluctuations in the cold dark matter (CDM) dominated structure formation scenario. The weak lensing 
distortion effect on the shapes of background galaxies due to the halos is a powerful means of measuring the mass and 
the mass profile independent of the dynamical state and the uncertain relation between baryonic and dark matter 
distributions (e.g. see [1] for a thorough review). Weak lensing is now recognized as one of the most promising 
methods for constraining cosmology including the nature of dark energy (e.g. [2]), which is the primary motivation 
for both on-going and planned high precision weak lensing surveys, such as the CFHT Legacy Survey [35], the Hyper 
Suprime-Cam Weak Lensing Survey [3], the Dark Energy Survey (DES) [36], and ultimately the Large Synoptic 
Survey Telescope (LSST) [37], Euclid [4], and the Joint Dark Energy Mission (JDEM) [38]. 

There is a promising weak lensing technique, the so-called stacked lensing or halo-galaxy lensing correlation, which 
has been increasingly recognized as a statistically robust method since the first measurement a decade ago ([5, 6]; 
also see [7-12] for the recent measurements). The halo-galaxy lensing correlation can be measured first by measuring 
the tangential cllipticity/shear component of background galaxies with respect to the halo center in each halo region 
(galaxy or galaxy cluster selected beforehand) and then by averaging the tangential shear profiles over all the sampled 
halos that are tagged by a similar halo richness indicator and are in the same redshift range. Although lensing 
information for individual halos is lost, the halo-galaxy lensing correlation has several notable advantages. First, the 
average lensing profile becomes insensitive to substructures and asphericity of the individual halos and also to the 
projection effect, i.e. the contamination arising from uncorrelated large-scale structure along the line-of-sight. This is 
because these "contaminating signals" are averaged out via the stacking, under the assumption that the Universe is 
statistically homogeneous and isotropic. The stacked shear profile thus allows to extract the halo mass enclosed within 
a three-dimensional sphere of a given radius, which is an important quantity when using the halo mass function for 
constraining cosmology. Second, since the signal probed is linear in an estimator of lensing shear or galaxy cllipticity, 
the halo-galaxy lensing is more robust to systematic errors: the systematic errors, if relatively well behaved, are 
to some extend averaged out, in contrast to the cosmic shear measurements, where the measurement rests on the 
two-point correlation based methods and therefore the systematic errors may accumulate. For example, if a relatively 
clean selection of background galaxies is available, which is usually the case for cluster lensing [12], the stacked lensing 
is not contaminated by the intrinsic alignment of galaxy images, which is one of the potentially dominant systematic 
errors in the cosmic shear measurement. 

Thus the halo-galaxy lensing allows a precise measurement of the mean halo mass and the mean mass profile, which 
is useful in calibrating the halo mass-observable relations [13]. The stacked shear profile can boost the signal-to-noise 
ratios at very small angular scales as well as very large scales, while only small solid angles or too small signals are 
available for individual halo lensing. The stacked lensing signals on small scales arise from the mass distribution 
gravitationally bound within each halo (the so-called one-halo term), while the large-scale signals arise from the 
surrounding mass distribution, such as filamentary structures, correlated with the sampled halos (the two-halo term) 
(see [7, 14] for the halo model based modeling). In particular, [15] recently studied that, for planned wide-area lensing 
surveys, the halo-mass correlation function at large angles can be in principle used to measure features of the baryonic 



acoustic oscillations and to probe the effect of primordial non-Gaussianity [16], which appear at projected radii greater 
than 100 Mpc corresponding to 10 degree scales for lensing halos at z ~ 0.2. 

In this paper we reformulate the halo-galaxy lensing correlation function in harmonic space, i.e. by employing the 
full-sky approach, since all the previous observational and theoretical works of halo-galaxy lensing have been done 
based on the flat-sky approximation (sec [17, 18] for the study of the lensing power spectrum based on the full-sky 
approach). We study how the full-sky correction is relevant on angular scales probed by future wide-area surveys, and 
is compared with the BAO features and with the effect of primordial non-Gaussianity. We indeed find a non-trivial 
scale-dependent correction to the flat-sky approximation, and that counter-intuitively, the correction is not confined 
to large angles. 

The structure of this paper is as follows. In Sec. II we reformulate the average tangential shear profile based on 
the full-sky approach. In Sec. Ill we compare the shear profiles computed from the full-sky and flat-sky approaches. 
Sec. IV is devoted to brief conclusions. Readers who want to skip the derivations could go straight to our main 
analytic results Eqs. (18) and (29) and then to Sec. Ill where we show our numerical results. Unless explicitly stated 
we assume cosmological parameters to match a flat, A-dominated CDM model that is consistent with the recent 
WMAP results [19]. 

II. FORMULATION 

A. Weak Lensing Shear on the Sky 

The bending of light by foreground mass, gravitational lensing, is described by a mapping of the position of a point 
on the sky with coordinates (a; 1 ,^ 2 ) (in the absence of lensing) to a different position with coordinates (x l ,x 2 ) (e.g. 
see [1] for a thorough review). We call x % the lens- or image-plane coordinates and i 1 the source-plane coordinates. 
The deflection or displacement vector on the sky is given by the gradient of a lensing potential ip: 

x* = x i - V, (1) 

where tp is a line of sight projection of the gravitational potential and ";" denotes a contravariant derivative with 
respect to the image coordinates. For a finite-size source such as a galaxy, lensing causes its image to be distorted. 
The distortion is described by differences in deflection angles between different parts of a source galaxy image. The 
small displacement vector connecting the two points in the observed (or lens-plane) coordinates, 8x l , is mapped from 
the vector in the source plane, 8x l , as 

S& = fa - ilr'.J) 6x> . (2) 

Although the above expressions are covariant and valid in any coordinate system, to define the shear components 
7i and 72, it is convenient to work in a locally orthonormal system, i.e. gij = Sij, where Sy is the Kronecker delta 
function. In such a coordinate system, where one does not have to distinguish between lower and upper indices, we 
define the shear by 
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The convergence is a particularly informative quantity because it is a direct line of sight integral of the matter density 
fluctuation field S m (see appendix A). 

The tangential shear around a halo is in the flat sky approach defined in a polar coordinate system with the halo 
at the origin. On the spherical sky, this needs to be generalized to a spherical coordinate system. Without loss of 
generality the halo center can be taken at the north pole due to the statistical isotropy of the Universe. In spherical 
coordinates, the metric is given by 

,2 



ds 2 = d6 2 + sin 2 9d4> 2 . (5) 



At each point on the sphere, this coordinate system defines a local basis (e^e^) which is orthogonal, but not 
orthonormal. Orthonormality can be realized by locally making the transformation to the coordinates (9,(p), where 
ip is defined by dtp = sm.9d<fi. We then define the tangential and transverse shears by in the (0, ip) system setting 

It = -7i) 7x = 72 (in (6,<p) system). (6) 

It will be useful to have expressions for k, jt and 7 X in terms of derivatives of the lensing potential with respect 
to 9 and <p. To derive these, we will first construct xp-ij in the (6, <fi) system and then use the fact that tp-ij is a rank 
two tensor to transform to the (6, ip) system. We thus start by calculating the covariant derivatives of i/j with respect 
to 9 and </>. The first covariant derivative is a derivative of a scalar and therefore simply equals the partial derivative: 
ip-.i = ip,i- The second derivative is the derivative of the covariant vector ip.^ and is given by 

1p;ij = 1p,ij - ^ij1p,k, (7) 

where the Christoffcl symbols are defined by 

Kc = 2-9 Qd (9dc,b + 9bd,c - 9bc,d) ■ (8) 

The non-vanishing Christoffcl symbols are 



T 9 ^ = -sin 9 cos 9 
„a „a cos 9 



(9) 



"* 1" sin9 
and all others Christoffcl symbols are zero. This gives 
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We now transform to the orthonormal system, where the basis is transformed as e^ • e^ ^ 1 — > e v • e v = 1, giving 

Using Eqs. (4) and (6), we can now easily identify 
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Eqs. (12) will be the starting point of the derivations of the relation between tangential shear and projected mass and 
the relation between tangential shear and the angular power spectrum of the halo-mass cross correlation. 

B. Tangential Shear 

A quantity of great cosmological interest is the angle- averaged tangential shear (7t)(0), where in general the angle 
average of a quantity a is defined as 

I C27T 



(a)(9) = —j d<Pa(e,4>). (13) 

It follows from rotational symmetry that after stacking, there is no loss of information if we restrict ourselves to the 
angle-averaged shear and convergence ("fr)(d) 1 (lx)(d) and (k)(9). 



A useful property of (7x)(#) and (^t){9) is that the former has to equal zero and the latter can conveniently be 
expressed in terms of the projected matter density. In the flat sky approximation [20-22], 

( 7T >(0) = «(< ff) - (k)(0), (14) 

where 

K{<e) = ^J de'e'{K){e') (15) 

is the mean convergence within a circular disk with angular radius 9. The convergence in turn is a direct line of sight 
integral of the matter overdensity (see appendix A) . 

We will now derive the full-sky equivalent of Eq. (14). Using Eqs. (12) and Eq. (13), we find 

2 \ smfl 

(«>(*) = ^ (W> + |v^W>)- (16) 

To relate the shear to the convergence, we need to define the mean convergence averaged within circular aperture of 
radius 9 (the full-sky version of Eq. (15)), 

i j-B p2-n 

"I JOl 



R(<6) = — ; / sin 9'd9 

V ; 2tt(1- cos 9)J 

1 f 

sm.9' d9'{n){9') 



1 — COSiy ./q 



2(1 -cos( 
Eqs. (16) and (17) above now allow us to write 
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This is the main result of this subsection. Note that in the small-angle limit, 9 <C 1, the flat-sky result (14) is 
recovered. 

C. Average Tangential Shear Profile: Halo-Mass Correlation Function 

An even more useful observable is the so-called " stacked tangential shear" or "halo-galaxy correlation function" , 
which is obtained by computing the tangential shear profile of background galaxy images relative to each halo [39] , 
obscrvationally represented by a galaxy or cluster, and then taking the average over all the sampled halos [5, 6, 23]. 
The stacked shear profile can be estimated as 

<( 7T »(0) = -L£( 7T )(0), (19) 

halos 

where Nh is the total number of halos averaged over. The ensemble average of this quantity gives the halo-shear 
correlation function or the halo-mass correlation function: 

«7t))(0) = (Sh(e')jT(O' + 0)) = (MO) (7t)(0)>, (20) 

where Sh(0) is the angular density fluctuation field of halos on the sky, and jt is defined relative to the position 
(9', </>'). The independence of 9' and <f>' follows from rotational symmetry, which comes from the statistical isotropy 
of the universe. 

By comparing the observed stacked shear to the value predicted by theory, one can learn about large scale structure 
and the background expansion history (e.g. see [7] for such a study). To calculate the theoretical ((7t)), one first uses 



the halo- mass cross power spectrum P hm (k, z) and the background cosmology to construct the angular convergence- 
halo spectrum Cf h (see appendix A). Then, in the flat sky approximation, this angular spectrum is used to calculate 
((jt)) as a function of 8 through 

«7T»(0) = J l £c?«MW), (21) 

where C ; is the angular power spectrum of cross-correlation between halos and the convergence field due to mass 
distribution surrounding the halos for a given source galaxy population, and J n is the n-th order Bcsscl function of 
the first kind. However, this relation is modified when the curvature of the sky is taken into account. Below we will 
derive the exact relation. 

We start by expanding the halo density fluctuation field and the lensing potential field in spherical harmonics: 

s h {e,ct>) = $>L*WM), 
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W,4>) = Y.at m Yim(e,<j>). (22) 
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The spherical harmonics are 



Y lm (0,4>) = f-^J^/rico^e^, (23) 

with PJ" the associated Legendre polynomials. The angular power spectra are in general defined by 

((a? m y4 m ,)=CrSu>6 mm , (24) 

The tangential shear can according to Eqs. (12) and (22) be expanded as 

IT = -\Y, ^m (YlmMO, <t>) ~ Y lmiVV ($, </>)) . (25) 

lm 

It is now straightforward to work out the stacked tangential shear: 

(( 7 t»(0) = <<5» 7 t(M)> 

= ~\ E Yf m (0) (Y lm . d8 (0, 0) - Y lm , vv {6, 0)) Cf h 
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lT,y^r( Y '"M0,cb)-Y m ., vv (9^)) Cf h (26) 



(using Yi m (Q) = 5 m .o y/(2l + l)/47r). The second (covariant) derivatives of a scalar in the (0,ip) coordinate system 
can be read off from Eq. (11) or (12), giving for the terms in parentheses 
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where Pi is the Legendre polynomial and P 2 is the associated Legendre polynomial with m = 2. Since k = -^y 2 "^ and 
V 2 Y lm = -1(1 + l)Y lm , 

Cr—j^C?. (28) 



Inserting Eq. (27) and Eqs. (28) into (26) gives the final result 

{(7T))(g) = E g r g ,_ 2 ^ l 1 1 ^ f(cc 
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In the small-angle, large-Z limit [17], !>1,S<1, 
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so that for small angles, the flat-sky result (21) is recovered. 

Eq. (29) is the main result of this paper as it gives the exact expression for the observed stacked shear profile in 
terms of the theoretical angular power spectrum of the mass-halo correlation. 



III. RESULTS: COMPARING FULL- AND FLAT-SKY APPROACHES 
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FIG. 1: Left panel: The average tangential shear profile (divided by bias) as a function of halo centric radii for lensing 
halos at redshifts Zh = 0.3 and 0.8. The solid curves are the results computed from the full-sky approach, developed in this 
paper, while the dashed curves show the results from the flat-sky approximation. The curves denoted by "Limber" show the 
results using the Limber approximation for the calculation of the halo-mass power spectrum (see text for the details). Right 
panel: Relative differences between various approximations and the exact full-sky result. Note that the difference is defined as 

<5{(7t))/<(7T>) = ( {lT )) flat /{{IT)) full - 1. 



Now that we have an exact, full-sky expression for the stacked shear profile, ((7t))(0) (see Eq. 29), we compare 
it to the result computed from the conventionally used flat-sky approximation. We do this by studying (("fr)) for 
halos at two lensing redshifts z% = 0.3 and 0.8. For the redshift distribution of source galaxies we assumed z ss 0.37 
for the functional form n(z) ex z 2 exp[— z/zq] given by Eq. (7) in [24] such that the mean source redshift (z) = 1. 
Then we used the galaxies sitting in the redshift range z = [1, 1.5] as for source galaxies, yielding the mean redshift 
(z) ~ 1.2. For the matter power spectrum needed to compute the angular power spectrum C; we used the non-linear 
power spectrum, obtained by applying the HaloFit prescription [25] to the linear power spectrum, for a concordance 
ACDM model that is consistent with the WMAP 7-year results [19]. For the halo bias we simply assume a linear 
bias throughout this paper: the halo-mass correlation strengths are different from the mass correlation strengths by 
some constant factor. Although a more precise modeling of the halo-galaxy lensing requires to use a suite of N-body 
simulations, our approach is adequate enough, because the primary purpose of this paper is to study the impact of 
full-sky approach being compared with the conventionally used flat-sky results. 



In Fig. 1 we study the stacked shear profile as a function of comoving transverse distances from the halo center, 
R = 9\h- Note that the comoving radial distances to halos at z/, = 0.3 and 0.8 are Xh = 840 and 1984 h~ 1 Mpc, 
respectively, for a flat ACDM model and therefore the comoving transverse length of 100 /i _1 Mpc corresponds to 6.82 
and 2.89 degrees, respectively, for the two lensing redshifts. We also study the effect of the Limber approximation [26], 
which is often used in the literature to calculate the angular spectrum of halo-mass correlations C\ lK (see Appendix A 
for the details). 

The left panel of Fig. 1 compares the average tangential shear profile, computed from the full-sky approach, 
with those computed from various approximations, while the right panel shows the relative differences. The three 
approximations we consider are: flat-sky shear with exact C^, flat-sky shear with Limber approximated C/"*, and 
full-sky shear with Limber approximated Cf K . The main interest of this paper is the comparison of the flat-sky 
approximation to the exact expression, so we focus now on the flat-sky with exact Cj case. The first thing to notice 
is that the error caused by the fiat-sky approximation is larger for Zh = 0.3 than Zh = 0.8, which can be explained by 
noting that a fixed transverse comoving distance to the line-of-sight corresponds to a larger angle for lower lensing 
redshifts. As for the scale-dependence of the difference, as expected, it starts at zero for R — 0. Then, it starts 
to deviate linearly from zero even at small distances, with the flat-sky expression underestimating the signal. The 
perhaps counter-intuitive errors at small angles arise from the leading order correction in the asymptotic approximation 
relating the 2nd-order Legendre polynomial Pf to the 2nd-order Bessel function in Eq. (30), as explained in detail in 
Appendix B. The difference reaches the 1% level around the BAO scale and then becomes greater at larger distances. 
The flat-sky expression overestimates the signal by more than 10% at R £ 400 /i _1 Mpc for Zh = 0.3. Thus we need 
to account for the full-sky effect for a %-level high-precision measurement of the lensing-halo correlation, expected 
from upcoming wide-area lensing surveys. 
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FIG. 2: The BAO feature in the stacked tangential shear profile (divided by bias) for zh — 0.3. We show the shear rescaled by 
a factor R to make the feature come out more clearly. The red and blue lines correspond to a variation of ±1% in the total 
matter density with the baryon density kept fixed. The magenta curve corresponds to u>b — 0, giving no acoustic feature at all. 
The dashed line is for the fiducial model, but using the flat-sky approximation. 



Recently [15] studied that the halo-galaxy lensing profile can be used to probe the BAO feature, appearing at 
R w 147Mpc, and to explore the effect of primordial non-Gaussianity, which become manifest at even larger scales. 
Since the error caused by the flat-sky approximation becomes significant at larger angles, it is useful to compare this 
error to these signals. 

In Fig. 2, we show the stacked tangential shear profile (divided by bias) for z\ % = 0.3, zooming in on the BAO 
feature. The figure shows that the flat-sky error is comparable to a change in signal due to a variation of about 1% 
in the total matter density uj m (i.e. 1% of the fiducial value for w m ), with the baryon density Ub being fixed. 

The local-type primordial non-Gaussianity parametrized by /nl causes a scale-dependent bias in the halo-mass 
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FIG. 3: The stacked tangential shear profile as a function of /nl compared to the flat-sky result. We assume a linear bias 
6i = 2. Left panel: The profiles themselves. Right panel: The difference relative to the /nl = 0, exact result. The error due to 
the flat-sky approximation is comparable to the effect of / = 0(10). 



correlation [16, 27-33]: the halo-mass cross power spectrum can be formally given as 

P hm (k) = [bt + Ab{k)} P mm {k) 1 
where P mm (k) is the matter power spectrum. The scale-dependent bias is given by 
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(31) 



(32) 



(33) 



which is roughly equal to the Hubble scale. In the above, 5 C = 1.68 is the critical overdensity for spherical collapse 
and b\ is the bias in the absence of non-Gaussianity. The function g(a) is the linear growth factor [34]: it is well 
approximated by the fitting formula 
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(34) 



with 7 « 0.55 (in GR, with a cosmological constant). 

Since Cf K is a projection of P hm (k) (see Appendix A), it follows from Eq. (31) that /nl affects the low I part of 
the angular spectrum and therefore the stacked tangential shear at large separations. Assuming b\ = 2, in Fig. 3 (left 
panel), we show the signal for /nl = 0, ±10, ±100 and compare it to the flat-sky signal for /nl = 0. In the right 
panel, we show the difference of non-zero /nl and of the flat-sky expression for /nl = 0, relative to the exact result 
for /nl = 0. The error due to the flat-sky expression corresponds roughly to a change in /nl of order 10. 



IV. CONCLUSIONS 



In this paper, we have developed a full-sky formalism for the halo-galaxy lensing correlation or the stacked lensing 
shear profile. We have updated the the flat-sky formula relating the tangential shear to convergence and found the 



exact expression for the stacked tangential shear (("/t}) in terms of the angular halo-mass cross power spectrum. For 
the stacked tangential shear, we have confirmed by direct calculation that errors caused by the flat-sky approximation, 
which is currently used in the literature, are very small on scales relevant for current data (< 1% for comoving distance 
R < 20h~ 1 Mpc). However, if in the future halo-galaxy lensing is used to constrain primordial non-Gaussianity by 
measuring the stacked tangential shear at separations of several hundreds of Mpc, the error due to the fiat-sky 
approximation becomes relevant, it being of order 10% around R = 400ft. _1 Mpc. 

The size of the errors due to the flat-sky approximation, as shown in Figs. 1, 2 and 3, can be compared with Fig. 6 in 
[15], which shows that the fractional errors in measuring the halo-galaxy correlation function at these large scales for 
an all-sky lensing survey are at about 20-30% level. Therefore the errors are well within the statistical measurement 
errors for existing and future lensing surveys. Nevertheless the exact, full-sky calculation is just as easy to do as the 
flat-sky approximation, our results can give a way to correct for or estimate the full-sky effect for future wide-area 
lensing surveys. 
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Appendix A: angular spectra and Limber approximation 
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FIG. 4: Left panel: The full-sky result for the angular power spectrum of convergence and halo cross-correlation (divided by 
the linear halo bias) is compared with the result obtained from the Limber approximation for lensing redshifts zh = 0.3 and 
zh = 0.8 (see Eqs. [A6] and [A7]). Right: The relative differences. 

To calculate {{"/t)), we need the angular cross spectrum of weak lensing convergence k with the (projected) halo 
density fluctuation field 5h- These quantities can both be expressed as line-of-sight projections of three dimensional 
quantities (the matter overdensity S m and the halo overdensity 5^ respectively): 



n(h) 



dxW K (x)Sm(xn,T Q - X ), 



(Al) 



and 



Sh(n) 



dxW h (x)S h (xn,T a -x), 



(A2) 
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with weight functions 

3 W 2o ..„ , ^ f°^ X*~X (n a {z„) 



W K ( X ) = ^m "™,o X (1 + *) / ^ ^-^ -^ , (A3) 

2 Jz Xs \ n a ,tot J 

and 

V^( x ) = ff(z)^, (A4) 

rih,tot 

where n s (z s ) is the distribution of source galaxies, J dz s n s (z s ) = n s< totj and nh{z) is the distribution of lens halos, 
J dzrih(z) = Uh,tot- Here \ is the comoving distance along the line of sight and we assume spatial flatness. 
After expanding these fields in terms of spherical harmonics, the exact angular cross spectrum is given by 

C? K = ^[dxW\x)[dx'W K (x')[d 3 kP hm (k,T -x,T -x')Mkx)Mkx% (A5) 

where P hm (k, r, r') 6(k' — k) = (5^(k, r) S m (k', r')}, and ji(x) denotes the l-th order spherical Bessel function. 

When we can obtain a cross-correlation of the lensing convergence with halos at a single redshift Zh (and comoving 
distance Xh), i-e. Tih{z)/nh,tot = 5{z — Zh), the above expression reduces to 

C/ 1 * = 2^2 / (Pkjxikxh) J d X ji(kx)W K (x)P hm (k,T - X h,r - x ). (A6) 

Expression (A6) is the one we use for our calculation of the exact halo-galaxy lensing correlation. 
In the Limber approximation, the expression simplifies to 



c hK . . H 2 WK {Xh) phm f l + 2 



TO - Xh 



Xh \ Xh 

3 /I + - 

2 



tf 2 O m , (1 + %) Xh 1 P hm (^ ' To - Xfc) (A7) 

We show both the exact angular and the Limber approximated spectra (divided by bias) for Zh = 0.3 and z% = 0.8 
in the left panel of Fig. 4. In the right panel, we show the relative difference of the Limber spectra with the exact 
spectra. For low I, the Limber approximation underestimates the spectrum by up to 12% (for Zh = 0.3). The difference 
becomes negligible at I ~ 25. 

Appendix B: Leading order correction for the flat-sky tangential shear profile 

The exact stacked tangential shear is given by 



and the flat-sky approximation gives 



<<7t))(*)= E^^-^^Pf (cos*) (Bl) 



l^c^Mie). (B2) 



The special functions involved are related by 

P/2(C ° S 6) = frff r% hm (1 + ^' 9)) ' (B3) 

where e is small and goes to zero as 9 —> 0, I — » oo. For small 9, the exact shear can thus be written 

«7T»W = ll^^-^P^osO) (B4) 

S^ < g + i)g-a)! , "' AWq } (1 + 0(g)) (B5) 

£ ^ ' -W) C, (l + ^r 1 + 0(r 2 )) | (1 + 0(e)) (B6) 

/dZ 7 rm\n f i , 3 /-i , nn-2- 



^j^wjc, M + -r i + o(r 2 ))^(i + o( e )). (B7) 
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Since the main contributions to the sum/integral come from I ~ _1 , the correction to the flat-sky expression of 
order l^ 1 inside the integral translates to a shear correction proportional to 9. This effect alone would imply that the 
flat-sky approximation underestimates the shear at very low 9. However, we also need to consider the 0(e) term. We 
checked numerically that for fixed I, e oc 9 2 at least as long as 16 <C 1. However, the main contribution comes from 
/ 9 ~ 1. We checked that for fixed I 9 = x of order unity, e = e(x/0, 9) oc 9 and e < so that this leads to an additional 
0(9) correction to the flat-sky shear. This effect on its own would imply that the flat-sky expression overestimates 
the true shear. We find numerically that the first effect wins (see Fig. 1) so that for low 9, the relative difference 
between the flat-sky result and the exact result is negative and evolves linearly in 9. 
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